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ABSTRACT 

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of 
182 stellar targets along the line of sight to the newly-discovered ‘ultrafaint’ object Reticulum 2 (Ret 2). For 
38 of these targets, the spectra are sufficient to provide simultaneous estimates of line-of-sight velocity (vios, 
median random error = 1.3 km s“’), effective temperature (Teff, = 464 K), surface gravity (logg, ^logg = 
0.54 dex) and iron abundance ([Fe/H], 5[Fe/H] = 0-45 dex). We use these results to conhrm 18 stars as members 
of Ret 2. From the member sample we estimate a velocity dispersion of = 3.6;!;Q g km s“' about a mean of 
(I'los) = 64.8;!;} 0 km s“’ in the solar rest frame (~ -90.9 km s“' in the Galactic rest frame), and a metallicity 
dispersion of cr[Fe/H] = 0.50 ;!;}}'[3 dex about a mean of ([Fe/H]) = -2.67 /q 34. These estimates marginalize 
over possible velocity and metallicity gradients, which are consistent with zero. Our results place Ret 2 on 
chemodynamical scaling relations followed by the Milky Way’s dwarf-galactic satellites. Under assumptions 
of dynamic equilibrium and negligible contamination from binary stars—^both of which must be checked with 
deeper imaging and repeat spectroscopic observations—the estimated velocity dispersion suggests a dynamical 
mass of M(/?h) ~ 5RhO’^,„y(2G) = 2.4/Q g x 10^ Mq enclosed within projected halflight radius /?h ~ 32 pc, with 
mass-to-light ratio « 2M(Rh)/Ly = 462/j57 in solar units. 

Subject headings: galaxies; dwarf — galaxies: individual (Reticulum 2) — (galaxies:) Local Group — galax¬ 
ies; kinematics and dynamics — methods: data analysis — techniques: spectroscopic 


1. INTRODUCTION 

The census of Local Group galaxies has been revised dra¬ 
matically and repeatedly over the last decade. Mining of the 
SDSS stellar catalog has yielded discoveries of ~ 15 low- 
luminosity, dwarf-galactic satellites of the Milky Way (e.g., 
Willman et al.|[2005j [Zucker et al.||20d^ [Belokurov et al.| 
2007|l. The PanDAS and PanStarrs surveys have found nearly 
two dozen new sate llites of M31 ( jMcConnachie et d)]|2009[ 
[Martin et al.||2013| ). All told, the pop ulation of known Lo ¬ 
cal Group galaxies has nearly tripled ([McConnachie[[20 1 2|l . 
Of the new members, perhaps the most intriguing are the 
Milky Way’s ‘ultrafaint’ satellites. These objects have low¬ 
ered the floor of the observed galaxy luminosity function from 
My -8 to My ^ -2, such that some galaxies’ t otal luminosi¬ 
ties ca n be dominated by a single red giant star ( [Martin et al.[ 
[2008)1. Moreover, the structural parameters of ultrafaints have 
blurred what was once an obvious distinction between the 
Milky Way’s dwarf-galactic satellites and its globular clus¬ 
ters. As a result, the proper classification of most ultrafaint 
objects now requires spectroscopic measurements of velocity 
dispersions, metallicities and metallicity dispersions that can 
indicate the presence of a dark matter halo. 

Most recently, the Dark Energy Survey (DFS) ha s revealed 
nine new Galactic satellites at southern latitudes (|Koposov[ 
[eTaLjllOBl [The DFS Collaboration et al.[[M5l ‘K15’ and 
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‘DFS 15’ hereafter). Seven of the new objects have sizes and 
luminosities characteristic of ultrafaints. One of them. Retic¬ 
ulum 2 (Ret 2), has already attracted attention for several rea¬ 
sons. First, Ret 2 is the nearest {D ^ 30 kpc) and most eas¬ 
ily detected of the newly-discovered objects. Second, Ret 2 
clearly has a flattened morphology, which may indicate ongo¬ 
ing tidal disruption or pe rhaps rotation. Third, using p ublic 
data from the Fermi-LAT,[Geringer-Sameth et al.[P015|l And 
evidence for gamma-ray emis sion that is consistent with dark 
matter annih ilation in Ret 2. [The Fermi-LAT Collaboration! 
[et al. ( [2015| l assign low significance to the gamma-ray signal 
based on unreleased Fermi-LAT data; however, in an indep en- 
dent analysis of the public data, [Hooper & Linden] ([201511 re¬ 
produce the original detection. In any case. Ret 2 is clearly an 
intriguing object, and the first question to settle is whether Ret 
2 presents chemo-dynamical evidence for dark matter. That 
is, is Ret 2 a globular cluster or a galaxy? 


Here we present results from an initial spectroscopic ‘re¬ 
connaissance’ of individual stellar targets along the line of 
sight to Ret 2. We identify a sample consisting of 18 member 
stars, which we use to characterize Ret 2’s chemodynamical 
properties. Specifically, we estimate the means and disper¬ 
sions of velocity and metallicity distributions, and we check 
for velocity and metallicity gradients that might provide clues 
about Ret 2’s dynamical state and formation history. Finally, 
we compare Ret 2’s properties to those of known dwarf galax¬ 
ies and globular clusters in order to determine which popula¬ 
tion can claim Ret 2 as its newest member. 


2. OBSERVATIONS AND DATA REDUCTION 

On 19 February 2015, we observ ed Ret 2 with the Michi- 
gan/Magellan Fiber System (M2FS, [Mateo et al.|2012[| on the 
6.5-m Magellan/Clay Telescope at Las Campanas Observa- 
tory, Chile. M2FS deploys up to 256 optical fibers—each with 
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an entrance aperture of diameter 1.2 arcsec—anywhere over 
a (f/11) focal surface with diameter 29 arcmin. The fibers 
feed twin spectrographs that can be configured independently 
in a wide variety of modes. For the Ret 2 observations, both 
M2FS spectrographs were configured identically in ‘HiRes’ 
mode, with order-isolation filters providing coverage over the 
range 5132-5186 A at effective resolution TZ ^ 18000. 

The left-hand panel of Figure shows a color-magnitude 
diagram for stars projected within! 5 arcmin of Ret 2’s center 
(<32000 = 53.9256 deg, S 2000 = -54.0492 deg), from the pho¬ 
tometric catalog derived by K15 in their analysis of public 
images from the Dark Energy Survey. With M2FS we tar¬ 
geted 182 stars (large circles in Figure[T]) selected as probable 
members of Ret 2 by the DBS’s Milky Way Science Work¬ 
ing Group, who also provided the coordinates necessary for 
designing M2FS plug plates (private communication, 2015). 
The right-hand panel of Figure shows positions of the tar¬ 
geted stars. We also allocated M fibers to regions of blank 
sky in order to estimate background. We observed the Ret 2 
field for a total of 2 hours (3 x 2400 s). At evening twilight 
the field was already rather low in the west, such that airmass 
during our observations increased from 1.3 to 1.7 and seeing 
hovered above 1.2 arcsec. 

We analyze the individual images (i.e. the single 40-min 
exposures) as well as the ‘stacked’ image consisting of the 
average over the three individual images. Other raw images 
include ancillary exposures of emission-line sources (‘arc’ 
spectra), continuum sources (‘quartz’ spectra), twilight expo¬ 
sures, bias and dark exposures, all obtained using the same 
configuration and detector binning (2 x 2) as the science ob¬ 
servations. 

We process the raw data using standard IRAF routines, fol¬ 
lowing procedures described in detail by Mateo et al. (2015). 
To summarize, we begin by performing overscan, bias and 
dark current corrections. We then average the three individ¬ 
ual exposures to produce a single, stacked science frame for 
each of the two spectrographs. From science frames corre¬ 
sponding to individual as well as stacked exposures, we ex¬ 
tract one-dimensional spectra. For a given spectrum we ob¬ 
tain the wavelength solution by fitting a 4*-order cubic spline 
to ^ 30 emission lines identified in the identically-extracted 
arc spectrum from the same aperture. Residuals to the wave¬ 
length solution typically have rms ^ 0.3 km s“'. We correct 
for variations in fiber throughput by dividing each spectrum 
by a normalized fit to the continuum in the twilight spectrum 
obtained in the same aperture. Finally, we estimate and sub¬ 
tract the mean contribution f rom sky background fo llowing 
the procedure introduced by |Koposov et al.| (pOTT), which 
avoids re-binning science spectra. 

In parallel, we also compute variance spectra that account 
for all known sources of noise in the processed science spec¬ 
tra. These include Poisson noise, readout noise, rms fluctua¬ 
tions in the bias and dark images, and background noise prop¬ 
agated from individual sky spectra into the mean sky spec¬ 
trum. We assign large variances (10'°°) to pixels affected by 
cosmic rays, which we identify as outliers above an iterative 
fit to the continuum in the science spectrum. 

Figure displays examples of sky-subtracted M2FS spec¬ 
tra (from stacked science frames) for nine representative Ret 2 
targets. Those in le ft-hand panels are from likely members of 
Ret 2 (Section |4.1[ ) and show the relatively weak absorption 
features that are typical of metal-poor subgiants. Spectra in 


upper right-hand panels belong to likely foreground contami¬ 
nants and show the stronger absorption characteristic of late- 
type dwarfs. Lower right-hand panels display spectra for two 
stars identified as blue horizontal branch candidates by both 
K15 and DES15. These are the two bluest stars in our sample, 
and while their continua are well sampled, they exhibit no ob¬ 
vious a bsorp tion features in the observed spectral range (see 
Section [33] l. Einally, Eigureshows two anomalous spectra 
that reveal their sources to be things other than singl e stars. 
We discuss these spectra in more detail in Section[331 


3. MODELLING OF M2FS SPECTRA 

We model the sky-subtracted M 2ES spectra following the 
procedure of|WaIker et al.|(|20151 ‘W15’ hereafter). Briefly, 
we adopt a spectral model of the form (see also Koleva et al.| 
|2009|[E)posov et al.|2011| i 


M(A) = P/(A)7’ A 


1-F 


Qm(^^ 3“ Flos 


( 1 ) 


where c is the speed of light, the polynomial Pi(X) gives shape 
to the continuum and 7’(A[l -F S'AM+uos j ^ ^ continuum- 

normalized template spectrum. The template is redshifted ac¬ 
cording to line-of-sight velocity, vios, and also by a second 
polynomial, Qm(X), that adjusts for systematic differences be¬ 
tween wavelength solutions of real and template spectra. We 
generate templates from a synthetic library that has previously 
been used for the SEGUE Stellar Parameter Pipeline (‘SSPP’, 
Lee et al.||2008a|b|l and which is calculated over a regular 
grid in effective temperature (Teff), surface gravity (logg) and 
metallicity ([Ee/FI]). 

Including nuisance parameters that specify polynomial co¬ 
efficients and smooth the templates according to the instru¬ 
mental line-spread function (LSF), the full model has 15 free 
parameters. Pour have physical meaning: vios, T^ff, logg and 
[Pe/H]. In order to obtain simultaneous estimates of all pa¬ 
rameters, we follow W15’s Bayesian analysis. We adopt the 
same likelihood function given by W15’s equation 9, and the 
same priors listed in W15’s Table 2, with one exception: for 
ho, the parameter that specifies the amount by which the tem¬ 
plate is smoothed in order to mimic the instrumental ESP, 
we adopt a prior that is uniform between 0.001 A and 0.1 
A (W15’s prior on is uniform between 0.05 A and 0.1 A). 
We find that this adjustment improves fits to the narrowest 
absorption features. 


3.1. Posterior PDFs 

We use the software packa ge MultiNest ( |Feroz & Hob-| 
|son||2008| |Feroz et al.||2009|l to scan the parameter space 
and generate random samplings of the 15-dimensional pos¬ 
terior PDF For a given parameter, X, we then calculate mo¬ 
ments of the marginalized, ID posterior PDF as follows. The 
first moment is the mean, X = which we take 

to be the central value. The_second moment is the vari¬ 
ance, a\ = (N— 1)“* 'Ylii=i(3i~X)'^, which we take to be the 
square of the Icr credibility interval. The third moment is 
skewness, S = N~^^-^^[(xi-X)/, which equals zero 
for a symmetric distribution. The fourth moment is kurto- 
sis, K = which distinguishes Gaus¬ 

sian distributions (K= 3) from ‘leptokurtic’ ones with sharper 
peaks and fatter tails (K > 3) and ‘platykurtic’ ones with 
broader peaks and weaker tails (K < 3). 
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Fig. 1. — Left: Color-magnitude diagram for stars within R <15 arcmin of Re t 2’s center jKoposov et al.|2015). The red line is the Dartmouth isochrone for 
age = 12 Gyr, [Fe/H]= —2.5, [Q/Fe]=-l-0.4, and m-M = 17.4 (Potter et aLpOOS) . Large circles enclose stars identified as probable Ret 2 members by the DES 
MW working group (private communication, 2015) and subsequently targeted by M2FS. Filled circles mark 38 stars whose spectroscopic measurements meet 
quality-control criteria. Of those, red circles identify stars whose spectroscopic properties are consistent with Ret 2 membership; black circles are spectroscopic 
non-members. Right: Standard coordinates for stars within 0.2 magnitudes of the isochrone shown in the left-hand panel. Markers and their colors have the same 
meanings as in the left-hand panel. The large shaded circle represents the M2FS field of view. 


For all 182 stacked M2FS spectra from Ret 2 targets. Fig¬ 
ure 1^ shows how moments of PDFs for physical parameters 
vary with median S/N per pixel. As S/N grows, the PDFs 
become narrower and more Gaussian (5 ~ 0, /T ~ 3). In sub¬ 
sequent analysis we consider only the sample consisting of 38 
observations for which the PDF for velocity has u < 5 km s“*, 
|5| < 1 and |A'-3| < 1 (red points in Figure]^. 

3.2. Accuracy and Precision 

In order to examine accuracy and precision, we apply the 
same analysis to 256 solar spectra that we acquired during 
evening twilight on the same night and with the same spec¬ 
trograph configuration used for the Ret 2 observations. Reas¬ 
suringly, parameter estimates from the collection of twilight 
spectra exhibit empirical scatter that is consistent with vari¬ 
ances of the posterior PDFs. Among the 256 twilight spectra, 
standard deviations of central values for vios, Fgff, logg and 
[Fe/FI] are 0.1 km s“’, 43 K, 0.04 dex and 0.02 dex, respec¬ 
tively. Furthermore, mean estimates show only mild system¬ 
atic offsets with respect to known solar values; (vios) - vios,© = 
-3.9 km s“', (Feff) - T^sq = 69 K, (logg) -logg© = 0.09 dex 
and ([Fe/H])-[Fe/H]Q = -0.20 dex. Following W15, for 
stellar-atmospheric parameters we treat these offsets as zero- 
point errors and subtract them from raw estimates obtained 
for science targets. We also add the empirical scatter, in 
quadrature, to (square roots of) second moments for all ob¬ 
servations. After applying these adjustments, our estimates of 
stellar-atmospheric parameters have median (minimum, max¬ 
imum) random errors of 177 ;^, = 464 (50, 1035) K, aiogg = 0.54 
(0.08, 1.30) dex and cr[Pe/Hj = 0-45 (0.04, 1.05) dex. 

For velocities, we do not adjust the zero point based on re¬ 
sults from the twilight spectra. In our experience, M2FS ve¬ 
locities shift systematically as ambient temperature changes 


rapidly during twilight, and the observed offset between twi¬ 
light and solar velocities reflects this phenomenon. As pointed 
out to us by the DES collaboration (private communication, 
2015), a temperature-dependent velocity shift continued, al¬ 
beit at a slower rate, during our observations of Ret 2. In order 
to examine this effect, we compare velocity estimates that we 
obtain from fits to spectra obtained in individual science expo¬ 
sures. We consider only the 73 single-exposure spectra whose 
results satisfy the same quali ty-control criteria imposed on the 
stacked spectra (Section [TT] i. Consistently with the effect no¬ 
ticed by the DES collaboration, we And that velocity estimates 
shift systematically over the three exposures and at different 
rates for the two spectograph channels. Erom 39 (34) spec¬ 
tra observed on M2ES’s ‘blue’ (‘red’) channel, we measure a 
weighted mean difference of (vexp 3 - Vexpi) = 2.1 ± 0.3km s“* 
(1.1 ±0.3 km s“') between the first and third exposures. Thus 
the raw velocities that we estimate from the stacked frames 
have channel-dependent zero-point errors. In order to com¬ 
pensate for this effect, we subtract the appropriate channel- 
specific value of (vexp 3 - Vexpi) from each of the raw veloci¬ 
ties that we estimate from the stacked frames. We also add 
the corresponding error in our estimate of (vexp 3 - Vexpi), in 
quadrature, to the error of our velocity estimates, after which 
our velocity sample has median (minimum, maximum) er¬ 
ror =1.3 (0.4, 3.3) km s“’. This procedure is designed 
only to remove the dependence of the velocity zero point on 
channel, and does not necessarily correct for a global error in 
zero point. Flowever, a comparison to an independent spectro¬ 
scopic sample obtained with VLT/ELAMES (Koposov et al., 
in preparation), shows that the zero points are in agreement. 
Eor the 17 stars common to both samples, the mean velocity 
difference is (vm 2 fs - vvlt) = 0.3 ± 0.3 km s“^. 
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Fig. 2.— Sky-subtracted M2FS spectra, with best-fitting models overplotted. Text gives target ID (upper left), equatorial coordinates (in degrees), color and 
magnitude, and estimates of redshift and stellar-atmospheric parameters. Panels on the left-hand side show spectra for probable members of Reticulum 2. The 
two upper-right panels show spectra for probable contaminants in the Galactic foreground. Lower-right panels show spectra for two blue horizontal branch 
candidates. 


3.3. Repeatability 

In order to gauge repeatability of our estimates for physi¬ 
cal parameters, we again consider results from our modelling 
of 73 single-exposure spectra whose results satisfy quality- 
control criteria. The top row of panels in Figure [^compares 
parameter estimates from the first exposure with mose in ei¬ 
ther of the two later exposures (all velocities from a given ex¬ 
posure and channel are shifted to remove the systematic drift 
in zero-point discussed above). The bottom row of panels 
in Figure displays distributions of deviations with repect 

to weighted means—(J) = i/^x)/Y!i2\ 

normalized by propagated error. Gray curves show Gaussian 
distributions with the same integrated area as the histograms. 
The observed distributions are all narrower than the Gaussian 


curves, implying that the variances of the PDFs give credibil¬ 
ity intervals that may be slightly overestimated. 


As a final check on reliability of our velocity estimates, we 
independently measure line-of-sight velocities using IRAF’s 
‘fxcor’ package. We cross-correlate each target spectrum 
from the stacked frame with a high-S/N M2FS spectrum that 
we acquired for the radi al velocity standa rd star CPD-432527 
(vios= 19.7±0.9kms“'; Udryetal. 1999 1 immediately before 


the Ret 2 observations. Reassuringly, we find good agreement 
with results from the Bayesian procedure described above. 
For the 38 spectra that satisfy quality-control criteria (Sec¬ 
tion 3.11, the mean deviation between (raw) velocity estimates 
is ^0.1 ±0.30 km s“'. Furthermore, when we model the 
spectrum of CPD-432527 using our Bayesian procedure, we 
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Fig. 3. — Sky-subtracted M2FS spectra for two anomalous sources. Top: 
This source is well-fit by the superposition (red) of two metal-rich model 
spectra (blue and green). Bottom: We speculate that the extremely wide, deep 
and irregular absorption features correspond to an absorption line system at 
redshift z > 2.5. 

obtain vios = 19.7 ±0.3 km s“', recovering the previously- 
published value. 


3.4. Data 

Table [T] lists results for the 38 spectra that satisfy the 
quality-control criteria discussed in Section [33 The online 
version of this article includes a supplementary table with re¬ 
sults for all 182 observed stars. The first three columns of 
Table [^list target ID and equatorial coordinates, followed by 
angular separation from Ret 2’s center and g, r magnitudes 
from the photo metry of K15, correcte d for extinction using 
the dust maps of|Schlegel et al.|(|1998|l. The seventh column 
lists median S/N per pixel. Columns 8-11 then list the re¬ 
sults of our spectroscopic modelling in terms of the first four 
moments of PDFs for vios, Teff, logg and [Fe/FI]. For each pa¬ 
rameter, the quoted central value is the first moment (mean) 
and the credible interval is the square root of the second mo¬ 
ment; both values have been adjusted acc ordin g to the offsets 
and empirical scatter described in Section [3721 Skewness and 
kurtosis are listed parenthetically as superscripts above the 
credible interval. The last column indicates whether or not 
the star satisfies membership criteria discussed in Section |47T] 
For all observations, the first exposure began at heliocentric 
Julian date HJD=2457072.542 days and the midpoint of the 
3 X 40-minute exposure occured at HJD=2457072.579 days. 


3.5. Special Cases 

Four spectra from our sample merit special attention. Table 
[^ lists their coordinates and g, r magnitudes. The first two 
(Ret2-161 and Ret2-034) belong to the two blue horizontal 
branch (BHB) candidates and are displayed in the lower-right 


panels of Figure [^ Their spectra show well-sampled con- 
tinua but no obvious absorption features, as expected given 
the greater degree of ionization at high temperature. As a re¬ 
sult, model parameters for these stars are poorly constrained. 
While their velociti es ar e loosely consistent with Ret 2 mem¬ 
bership (see Section |47T] l, the associated uncertainties are ~ 30 
km s“' and ~ 260 km s“', such that nei ther spectrum satis¬ 
fies our quality-control criteria (Section |3.1| l. On the other 
hand, their temperatures (7’eff=6890 ± 680 K and 6880 ±630 
K), while also highly uncertain, are the two largest that we ob¬ 
tained for any of our 182 observed spectra. We conclude that 
these two stars remain strong BHB candidates, with confirma¬ 
tion requiring spectra that cover a larger range in wavelength. 


The next special case (Ret2-097) appears in the upper panel 
of Figure This ‘star’ has color and magnitude placing it 
near the isochrone and above the horizontal branch in Fig¬ 
ure However, the spectrum is rich in broad and apparently 
double-valleyed absorption features. We model this spectrum 
as the superposition of two spectra of the form given by Equa¬ 
tion!^ Allowing for a second set of stellar-atmospheric pa¬ 
rameters and a second continuum polynomial, the resulting 
model has 25 free parameters. As shown in Figure [^ this 
double-star model gives a reasonably good fit to the spectrum. 
One of the sources has vios = 54.8 ±0.6 km s“\ slightly be¬ 
low the range of velocities that we attribute to Ret 2 members 
(Section [ 4 . Ij i. Its strong surface gravity (logg = 4.41 ±0.21) 
and rich metallicity ([Fe/H] = -0.37 ±0.08) are typical of 
foreground contamination. The other contributing star has 
Vios = 120.4 ±0.7 km s“\ logg = 4.52 ±0.23 and is also 
metal-rich, with [Fe/H] = -0.52 ±0.10. We speculate that 
this source is a physical binary composed of two K-type main 
sequence stars at a distance of a few kpc. Assuming a mass 
ratio near unity, the center of mass has vios ~ 100 km s“\ in¬ 
consistent with Ret 2 membership. Due to the anomalous na¬ 
ture of this spectrum, we do not include these results in Table 
[T]or consider them in subsequent analysis. 

Finally, the bottom panel of Figure displays a spectrum 
(Ret2-019) that is difficult to classify given the limited spec¬ 
tral range. Here we consider a few possibilities. First, the ab¬ 
sorption features may have molecular origin, specifically from 


MgH, which has a bandhead at ~ 5200 A and extends blue- 
ward by ~ 200 A. However, the features in our spectrum are 
too wide, too deep and too irr egularly spaced to m atch the typ¬ 
ical appearance of that band ( Week et al.|2003 1 . In addition, 
MgH absorption becomes evident typically in stars of spec¬ 
tral types KO or later, much redder than the de-reddened color 
ig-r ^ 0.17) of this source. Alternatively, this spectrum may 
reveal the presence of a complex absorption line system along 
the line of sight to a remote active galactic nucleus. The line 
widths and extreme de pths appear similar to numerous line 
complexes reported byjSrianand et al.|(|2010|l for various ab¬ 
sorption lines in Lyman-a systems. Furthermore, the source’s 
blue color is consistent with that of a background quasi-stellar 
object (QSO) at moderate redshift. For a number of plausible 
UV lines that may produce the observed absorption features, 
the implied redshift of the absorber is z > 2.5, with the back¬ 
ground source obviously more distant. The nearest known 
extragalactic source in the NED catalog is detected in the UV, 
by GALEX, and is 24 arcsec away from our target. A defini¬ 
tive classification of this source will require spectroscopy with 
broader wavelength coverage. 


4. RESULTS 





















6 


_ 17 
M 18 

a 19 
^ 20 
21 

CB 100 

• !—I 

03 

O 

t: 10 

ri 

CB 1 
03 C 
03 ^ 

S n 
^ 0 

^ -5 

(M C® 


:/ 


t’lc 


» V** • * 


q; 100 

S 10 

ce 1 

•j: 1 

$ 0.1 


I llllllll I lllllll Jill 


. 0 ' 


= - - 100 


5 = 




iii m iii llllllll iiiiiiii i 


r'V' ■ ■ w 


. % ^ 

S[km/s]*- 

I IIIIIIM I lllllll I III 




mm 


0 

-5 


t-h O O 
t-h O 



^ O 


rH O 


^ O 


S/N 


S/N 


S/N 


S/N 


Fig. 4.— Moments of marginalized, ID posterior PDFs of model parameters vs median S/N per pixel, for all 182 spectra acquired wit h M 2FS. Histograms 
show marginalized distributions for each moment as well as S/N. Red points indicate observations satisfying quality-control criteria (Section [3.1) , where posterior 
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Gaussian distributions (gray curves) with integrated areas equal to those of the histograms. The empirical distributions are generally narrower than the Gaussian 
ones, implying that our en'ors may be slightly overestimated. 
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TABLE 1 

M2FS Stellar Spectroscopy of Reticulum 2“ 


ID 

0:2000 

[hh:mm:ss] 

<^2000 

R 

[arcmin] 

g 

[mag] 

r 

[mag] 

S/N* 

Vlos 

[km s“^]^ 

7'eff 

[K] 

logio^ 

[dex]'' 

[Fe/H] 

[dex] 

member? 

Ret2-004 

03:35:56.27 

-54:03:16.2 

2.1 

18.83 

18.37 

5.7 

63.9±2.3''’-^’^''> 

6506±8iT^'3'^'^> 

3.12± 1.14‘-°-’''^''') 

-2.26±0.64<^A2.6) 

Y 

Ret2-006 

03:36:01.76 

-54:04:05.4 

3.1 

19.52 

19.05 

3.2 

62.5± 1.9*"’-^’> 

5313±524<“'‘’^'^) 

4.02±0.62L‘‘.4.9) 

-2.88 ± 0.50<“'^’^ "> 

Y 

Ret2-016 

03:36:35.74 

-54:01:56.9 

7.9 

17.99 

17.58 

9.5 

133.9 ±0.5*““’^“> 

5460±78<““'^‘’' 

4.56±0.14*“'’^'’' 

-0.82±0.06‘“““’^'’ 

N 

Ret2-021 

03:36:22.83 

-53:59:55.5 

6.7 

18.53 

18.07 

7.7 

128.3 ± 0.5<^"-^ ’'’ 

5285± 115<“’^-^* 

3.19±0.15<“'^-^" 

-0.99 ± 0.12*“^’^ 

N 

Ret2-023 

03:36:21.86 

-54:00:40.6 

6.3 

20.23 

19.81 

1.6 

66.7 ± 2.0<"'^-'‘“* 

5658 ± 1034<‘’'^'^"> 

3.37±0.89*^-^'^'‘’ 

-1.96± 1.05<““‘-2-‘’> 

Y 

Ret2-032 

03:35:40.69 

-54:10:05.1 

7.1 

18.19 

17.90 

4.9 

49.0 ± 1.0'^“‘‘> 

5786±503<“’^’* 

3.81 ±0.48<““’> 

-1.75±0.45<“’^'^' 

N 

Ret2-035 

03:35:37.06 

-54:04:01.2 

1.3 

18.54 

18.03 

4.6 

62.9± 1.2*"‘‘-^®> 

5034±441<“'’’'‘'’* 

1.50±0.89‘““-^> 

-2.90±0.45*“-^’^-^' 

Y 

Ret2-040 

03:35:58.14 

-54:02:04.7 

2.5 

19.25 

18.80 

4.1 

69.1±1.5<-“-^-^‘) 

6283 ± 865*““ ‘'^-^' 

2.97± 1.05*^‘‘-^''> 

-2.28±0.71<““’^-^> 

Y 

Ret2-049 

03:35:20.36 

-54:18:16.9 

15.7 

20.55 

20.26 

1.2 

308.6 ±2.3<“‘’-^’> 

6009±645<“'’^** 

4.06±0.59<-"‘’-^*> 

-0.58 ± 0.67<““‘1.3.i) 

N 

Ret2-062 

03:35:09.50 

-54:02:29.6 

4.8 

18.20 

17.73 

7.3 

95.9±0.6'“^‘-^"> 

5426 ±98'“''^** 

4.49 ± 0.17'“'^-^*> 

-0.87 ± 0.09*“ 

N 

Ret2-063 

03:35:11.65 

-54:03:21.8 

4.5 

19.28 

18.80 

3.0 

70.4± l.9*^-^’^-^> 

4494 ±548<^'’‘’-^* 

I.19±0.85<“-^'^> 

-3.52±0.52<''^’^’> 

Y 

Ret2-069 

03:35:02.49 

-54:03:54.2 

5.9 

19.20 

18.72 

3.7 

67.9± 

5048±523<“*’''‘* 

1.90 ± 0.83'“^ 

-2.77 ± 0.55<“'^’^'’' 

Y 

Ret2-072 

03:34:37.97 

-54:13:59.6 

14.5 

20.47 

20.06 

1.4 

137.3 ±3.0*“'^’^-'” 

5574±762<“’^‘> 

3.39 ± 1.20<^‘’-^“> 

-2.13±0.68<“'^’^’> 

N 

Ret2-076 

03:34:36.69 

-54:06:44.9 

10.3 

18.84 

18.30 

4.3 

71.6±0.7‘""-^‘” 

4988 ± 229<“ ’’^ ‘’* 

2.26 ± 0.29'“ ''^'^' 

-1.52±0.26<“’^-^’ 

Y 

Ret2-077 

03:34:53.23 

-54:14:03.9 

13.2 

20.63 

20.26 

1.2 

290.8 ± 2.7<^'‘-^ ‘‘' 

5372± 464<“’^‘* 

4.28 ± 0.50‘“‘ ‘‘-'^ *> 

-1.57±0.55'“““*’ 

N 

Ret2-079 

03:34:54.24 

-54:05:58.0 

7.6 

18.91 

18.42 

3.9 

70.0± 

4915±543<“’’^'^> 

3.44 ± o.74<^'^-^'‘> 

-3.19±0.52<“'^’^‘’> 

Y 

Ret2-080 

03:34:47.93 

-54:05:25.0 

8.3 

17.50 

16.91 

11.2 

63.2±0.5‘"‘-^-^* 

4501 ± 141<“’^'^> 

0.58 ± 0.48*“ ’-^ '^> 

-3.16±0.16'“’’^'^’ 

Y 

Ret2-087 

03:34:13.05 

-53:59:56.1 

13.4 

19.50 

19.00 

2.9 

200.5 ± 1.7<^'‘-^‘’> 

5211± 150’““’^“' 

4.55 ± 0.26'^ *-^ ‘'' 

-1.32±0.18*“'’^“' 

N 

Ret2-089 

03:34:05.48 

-54:03:49.9 

14.2 

17.61 

17.01 

10.1 

53.2±0.5*^"'^’> 

5107±50<-"‘‘-^'> 

4.77 ± 0.1 

-0.49 ± 0.04<“ 

N 

Ret2-090 

03:33:41.70 

-54:00:07.2 

17.9 

18.30 

18.21 

4.5 

32.4 ± 1.6'"'^’^-^> 

6465 ± ni^-o.2.23) 

3.43± 1.30‘“‘'^-^‘'' 

-1.74±0.62<““'='^’> 

N 

Ret2-092 

03:34:49.20 

-53:50:19.7 

14.8 

17.55 

16.98 

11.3 

16.2±0.4'""-^-’> 

5218±58'“-^'’* 

4.60±0.12<“'^-^*> 

-0.01 ± 0.05*“ “’^ 

N 

Ret2-100 

03:34:39.81 

-54:00:58.5 

9.4 

18.18 

17.95 

8.1 

199.4 ± 0.9*^ 

6368±333<“'’^® 

4.05 ± 0.36'^"-^ ''> 

-1.07±0.24'““-‘-^-*' 

N 

Ret2-102 

03:34:18.31 

-54:10:06.2 

14.2 

18.06 

17.47 

7.5 

28.1±0.5<‘^‘-^"> 

5118±58<^'‘-^“> 

4.70±0.13<^'‘-^'^> 

-0.54 ± 0.06<““ “-^ “> 

N 

Ret2-103 

03:34:13.93 

-54:09:34.4 

14.5 

18.33 

17.71 

6.3 

20.4±0.5‘"‘-^“’ 

5069±52<^'^'^-^> 

4.79 ± 0.10<““’-^ *> 

-0.47 ± 0.06‘““ “-^ 

N 

Ret2-110 

03:34:29.93 

-54:11:11.7 

13.4 

19.80 

19.35 

2.4 

290.5 ± 1.9<^‘-^*> 

5797 ± 888<“-^’^-^* 

2.12± 1.1 

-1.69±0.82<“'’^'‘> 

N 

Ret2-Ill 

03:34:30.06 

-54:09:22.2 

12.4 

18.20 

17.88 

8.0 

270.1 ±0.7<^‘-^‘’> 

5523± 182<“'^’^'‘* 

3.63±0.23<“'^'^'^> 

-1.22±0.17<“'^’^‘> 

N 

Ret2-115 

03:35:31.13 

-54:01:48.2 

2.0 

17.57 

17.01 

9.0 

60.4±0.7<^'^’^'^) 

5164±242<“’^'^) 

2.44 ± 0.63*“‘'’-'^ “> 

-2.62±0.26<“'’^‘’> 

Y 

Ret2-126 

03:35:20.96 

-54:03:48.1 

3.2 

18.91 

18.43 

5.0 

65.6± 

5185±572<“*’‘'‘> 

2.91 ±0.62*““’^'^> 

-2.67 ± 0.57<“’^ ‘’> 

Y 

Ret2-128 

03:35:06.56 

-54:06:04.3 

6.1 

19.74 

19.26 

2.3 

160.5 ± 1.6<“-^'^‘‘* 

5466±319*“-^’”' 

3.99 ± 0.54<^'^-^ ''> 

-1.02±0.33'“““'” 

N 

Ret2-129 

03:34:57.57 

-54:05:31.4 

7.0 

18.90 

18.40 

3.0 

61.8± 1.4<"'^'^-^> 

5296±539<“'^’^'^> 

2.99 ± 0.65*^ 

-2.04±0.57<““'’2-'’> 

Y 

Ret2-134 

03:35:15.17 

-54:08:43.0 

7.0 

19.71 

19.26 

3.0 

70.2±3.3*"“-^> 

5875± 1019‘"‘'^‘’> 

1.94 ± 1.15<“”'^> 

-2.41 ±0.87<““''^'> 

Y 

Ret2-136 

03:35:14.01 

-54:05:58.1 

5.1 

19.98 

19.55 

1.8 

66.4±2.9‘"-^-^'^' 

5582± 849<“‘‘’^'^* 

2.78± 1.25'^-^-^-^’ 

-2.14±0.84<““’^'‘' 

Y 

Ret2-138 

03:35:13.73 

-54:04:56.7 

4.6 

19.67 

19.22 

1.9 

60.1±2.1<"“*> 

5032±694<“’’^‘> 

2.97± 1.02<^‘1.2.7) 

-2.62±0.72<“’^'^> 

Y 

Ret2-142 

03:35:44.18 

-54:01:50.0 

1.2 

20.34 

19.89 

1.7 

65.0±2.9'^'^-^*> 

6067±529<““‘'^*> 

4.39±0.41<“‘'^-‘‘'^> 

-1.24±0.60‘““'^'^“> 

Y 

Ret2-147 

03:35:36.93 

-53:54:45.1 

8.2 

20.52 

20.15 

1.5 

302.3±2.1<^'^-^‘> 

5335 ± 524<“'^’^'^* 

3.71±0.68<^'‘-^‘’> 

-1.23±0.58*“'’^‘> 

N 

Ret2-153 

03:35:46.16 

-54:07:33.8 

4.6 

19.21 

18.74 

4.3 

217.9 ± 1.2*“-^’^'* 

5089 ± 333<“''’‘*“> 

2.82±0.5l‘^-^''‘-^’ 

-2.16±0.37*“'’^‘> 

N 

Ret2-159 

03:35:49.95 

-54:03:21.4 

1.2 

19.65 

19.22 

3.1 

293.2± l.l<“-^"> 

5329 ± 235<“'^’^ '’* 

4.26±0.37<^'‘-^'^> 

-1.37±0.25<“'’^'’' 

N 

Ret2-178 

03:36:07.75 

-54:02:35.5 

3.8 

17.39 

16.81 

15.6 

61.8±0.4<^‘’^") 

4688± 105<“'’^‘’* 

1.34±0.27<^*-'‘'^> 

-2.74±0.13*“““’^-*' 

Y 


"This version lists results only for observations satisfying quality-control criteria (Section [3.1| . The online version includes results for all targeted stars, 
^median signal-to-noise ratio per pixel 
“^line-of-sight velocity in the heliocentric rest frame 
'^units of g are cm s~' 
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TABLE 2 

Coordinates and photometry eor special cases (Section |3.5> 


ID 

Otiooo 

[hh:mm:ss] 

^2000 

R 

[arcmin] 

s 

[mag] 

r 

[mag] 

notes 

Ret2-161 

03:36:18.68 

-53:57:45.1 

7.5 

17.99 

18.17 

BHB canididate 

Ret2-034 

03:35:39.85 

-54:04:58.1 

2.0 

18.28 

18.58 

BHB canididate 

Ret2-097 

03:34:37.33 

-53:53:53.9 

13.1 

17.66 

17.10 

double star 

Ret2-019 

03:36:22.99 

-53:55:07.9 

9.9 

21.11 

20.94 

unknown 


Scatterplots in Figure show how the spectroscopically- 
derived quantities (vios, Feff, logg, [Fe/H]) relate to each other 
and to photometrically-derived quantities (position, color, 
magnitude). 


4.1. Membership 


Figure]^ also helps to distinguish bona fide members of Ret 
2 from contaminants in the Galactic foreground. We expect 
Ret 2 members to exhibit a relatively narrow range in veloc¬ 
ity, lower metallicities, lower surface gravities, and to be clus¬ 
tered nearer Ret 2’s center than are their foreground counter¬ 
parts. Figure 1^ clearly shows a population of stars having 
these characteristics. Red boxes in Figure [^enclose measure¬ 
ments that cluster near the centers of Ret Ts distributions for 
Vios, logg, [Fe/H] and position, as determined by eye. We 
make no selection based on color, magnitude or temperature 
because the full ranges for these quantities are consistent with 
Ret 2 membership. We count 18 stars that are inside all red 
boxes. Table[2lists a membership status of ‘Y’ for these stars, 
and ‘N’ for stars that lie outside any of the red boxes. 


In previous work we have employed an expectation- 
maximization (EM) algorithm in order to es timate model 
parame ters in the presence of contamination ( [Walker et al.| 
|2009b| ). Under the assumption that velocity dispersion is spa- 
tially uniform, our EM algorithm provides estimates of mem¬ 
bership probabilities for all individual stars. Reassuringly, 
when we apply the EM algorithm to the Ret 2 data in Table 
the sum of membership probabilities is 17.8, with each of the 
stars inside our red boxes receiving membership probability 
> 0.85. Eurthermore, the EM algorithm’s estimates of means 
and dispersions of velocity and metallicity distributions agree 
well with those obtained in the independent analysis that we 


present in Section 4.2 


4.2. Chemodynamics of Reticulum 2 

In order to characterize the internal chemodynamics of Ret 
2 , we consider a simple model in which the velocities and 
metallicities of the 18 member stars are normally distributed 
about means that can vary systematically with position in or¬ 
der to account for velocity and metallicity gradients. Specif¬ 
ically we assume that the mean line-of-sight velocity, is 
a function of projected position {R, 9), 

tiv,JR,9) = (vios) +k„,„Rcos(6»- (2) 

where (vios) is the mean velocity at the center, = 
|d/i,,|^^/dR| is the magnitude of maximum velocity gradient 
and (measured from north of Ret 2’s center and opening 
to the east, in equatorial coordinates) specifies its direction. 
We assume that mean metallicity, /i[Fe/H]> is a function only 
of separation from the center, 

M[Fe/H](-K) = ([Fe/H]) -|-k[Fe/H]R, (3) 


where ([Ee/H]) is the mean metallicity at the center 
and k[Fe/H] = dF[Fe/H]/d^ specifies the magnitude of any 
(isotropic) metallicity gradient. 

Under these assumptions, the joint probability density for 
observables vios and [Ee/H], for a star at position {R, 9), is 


p(vios,[Ee/H]|R,0) 


(27r)-' 


X exp 


1 (vios-/Tv,„)^ _ 1 ([Ee/H]-/r[Fe/H])' 
2 ' 


2 “I” 

^ ^ Vios ^ ^Vios 


2 r2 

[Fe/H] “'""[Fe/H] 


(4) 


where (Jv,„ and (T[Fe/H] are velocity and metallicity dis¬ 
persions, respectively, and and (5[Fe/H] are measure¬ 
ment errors. Given a vector of free parameters 9 = 
((t'ios),o-v,j^sJ^i’iosJ^''iosi ([Fe/H]),cr[Fe/H]i^[Fe/H])^ a data set con¬ 
sisting of N observations of Ret 2 members, D = 
{(R/,^,, Vios,;, [Ee/H],.)}(^p has likelihood 


C(D\9) = nXv'ios/, [Ee/H],.|R,-,0,). 


(5) 


Erom Bayes’ theorem, given the data, the model has posterior 
probability distribution function (PDE) 


p{9\D): 


CiD\9)pi9) 
p{D) ’ 


(6) 


64.8/|;i km s 
3.6/og km 


where p(9) is the prior PDE and p(D) = / C(D\9)p(9)d9 is 
the ‘evidence’. 

Again we use MultiNest to scan the 7-dimensional parame¬ 
ter space and to draw random samples from the posterior PDE. 
Eigurej^displays random samplings from PDEs for the means 
and dispersions. Table summarizes the results. Eor each 
of the seven free parameters, the second column lists bound¬ 
aries inside which the prior is uniform and nonzero (outside 
these boundaries, the prior probability is zero). The third col¬ 
umn lists median-likelihood values and credibility intervals 
that enclose the central 68% (95%) of area under the posterior 
PDE. 

We find that Ret 2 has mean line-of-sight velocity (vios) = 
with a resolved velocity dispersion of = 
s *. We obtain a mean metallicity of ([Ee/H]) = 
-2.61^W, with a dispersion of CT[Fe/H] = 0.50 /q} 3. Taken at 
face value, the estimated metallicity dispersion is five times 
larger than is exhibited by any known globular cluster (jWill-j 
[man & Strader||2012| ). However, we note that the estimated 
metallicity dispersion is also similar to the median credibility 
interval for our metallicity estimates. [Koposov et al.| ( p()ll| l 
demonstrate that, in such cases, estimates of intrinsic disper¬ 
sion can be extremely sensitive to systematic over- or under¬ 
estimation of measurement errors. Given the beh avior in the 
bottom-right panel of Eigure|^(see Section 3.3 for discus¬ 
sion), we suspect that our met^icity errors are systematically 
overestimated, which would imply that our estimate of Ret 
2’s intrinsic metallicity dispersion is underestimated. In any 
case, this estimate should be interpreted with caution and will 
need to be confirmed using spectra spanning a larger range of 
wavelengths. 

Einally, we do not detect significant gradients in either ve¬ 
locity or metallicity. Our estimate of = O.T/g^km s“' 
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Fig. 6.— Scatterplots showing relations among photometrically- (magnitude, r, color, g — r, and separation, /?, from Ret 2’s center) and spectroscopically- 
derived (vios, Fgff, logg and [Fe/H]) quantities for individual Ret 2 stars. Red boxes are drawn, by eye, to enclose probable members. For each observable, 
histograms show ID distributions for the full sample (open black) and likely members (solid red). 


TABLE 3 

Summary of probability distribution functions for chemodynamical parameters 


parameter 

prior 

posterior 

description 

(vios) [kms“‘] 

uniform between -500 and -1-500 


mean velocity at center 


uniform between 0 and -f500 

^ ^-K).9(+2.1) 

velocity dispersion 

[km arcmin”*] 

uniform between 0 and -FlO 

p, .-K).4(+1.0) 
^■‘*-0.2(-0.3) 

magnitude of maximum velocity gradient 


uniform between -180 and -1-180 

.,q+220(+251) 

'^-55(-93) 

direction of maximum velocity gradient 

([Fe/H]) 

uniform between -5 and -l-l 

- ^.y+0.34(+0.72) 
^•”'-0.34(-0.67) 

mean metallicity at center 

O'lFe/H) 

uniform between 0 and -f2 

U.JU_Q j3(_q 23) 

metallicity dispersion 

^[Fe/H] [dex arcmin“*] 

uniform between -1 and -l-l 

^ ^.^+0.06(+0.12) 
-0.06(-0.13) 

magnitude of metallicity gradient 
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Fig. 7.— Samples drawn randomly from posterior PDFs for means and 
dispersions of Ret 2’s metallicity (top) and velocity (bottom) distributions. 
Histograms display mai'ginalized, ID PDFs for each parameter. 

arcmin“* is consistent with zero at the < 2a level. On the 
other hand, values as large as ~ 1 km s“^ arcmin“* are sim¬ 
ilarly allowed. The steepest allowed gradients are directed 
along position angle ^ -80 degrees, roughly ~ 30 degrees 
from Ret 2’s morphological major axis (K15, DES15). Our 
estimate of A:[pe/H] = 0.02j;QQgdex arcmin“' is consistent with 
zero at the ^ Icr level, but also allows gradients as steep as 
±0.15 dex arcmin“' within ~ 2cr. We notice degeneracies be¬ 
tween our estimates of dispersions and gradients for both ve¬ 
locity and metallicity distributions, such that larger gradients 
correspond to smaller dispersions. However, our estimates 
for individual parameters (Table 1^ naturally include such ef¬ 
fects, as the ID PDFs for individual parameters are obtained 
by marginalizing over all other dimensions of the parameter 
space. 


4.3. Scaling Relations 


K15 use photometric data to show that Ret 2 occupies a 
region of structural parameter space that is populated by ob¬ 
jects of ambiguous classification, somewhat intermediate be- 
tween well-establish ed globular clusters and dwarf galaxies 
(|Gilmore et al.|2007| cf. K15’s Figure 17). Specifically, with 
projected halfiight radius = 32 ± 1 pc. Ret 2 is larger than 
nearly all globular clusters and smaller than nearly all known 
galaxies. Moreover, its absolute magnitude. My = -2.7 ± 0.1, 
would place Ret 2 among the least luminous members of ei¬ 
ther population. In these regards. Ret 2 and many of its newly- 
discovered siblings are similar to ‘ultrafaint’ satellites Segue 
1, Segue 2 and Willman 1, as well as to the extended globu- 
lar clusters Pal 14 and Crater (|Belokurov et al.|2014[|Laevens| 
let al.|2014| Mateo et al. 2015). 


Our spectroscopic results provide new leverage that can set¬ 
tle the question of Ret 2’s nature. For Galactic globular clus¬ 
ters as well as the dwarf spheroidal satellites of the Milky 
Way and M31, the top panel of Figure l^plots mean metallic¬ 
ity against luminosity. While globular musters show no obvi¬ 
ous trend, it is well-known that dwarf galaxies follow a lumi- 
nosity/metallicity relation ( |Mateo|[T998[ |Tolstoy et al.||2009{ 
|Kirby et al.||2013| l. Given the low mean metallicity that we 
estimate from the M2FS spectra, ([Fe/H]) = -2.67 /q 34, we 
place Ret 2 squarely onto the galactic relation (large black 
square in Figure]^. 


The lower panel of Figure ^ shows another well-known 
scaling relation that distinguishes dwarf galaxies from glob¬ 
ular clusters. Specifically, the mass-to-light ratios of dwarf 
galaxies are anti-correlated with luminosity, such that grav¬ 
itational potentials in the least l uminous galaxies all seem 
to be dominated by dark matter (|Mateo et al.|: 


|et al.|[2007l |Martin et al.||2007j Simon & Geha||2667| l. In 


1993 


contrast, the stellar kinematics of globular clusters generally 
do not require an internal dark matter component. In order 
to see these trends, we plot (GTy) against luminos¬ 

ity, where is halfiight radius, a^^^ is line-of-sight velocity 
dispersion, Ly is total V-band luminosity and G is Newton’s 
constant. The combination of macroscopic observables on 
the vertical axis has dimensions of mass-to-light ratic^ and 
is therefore sufficient to highlight the different behavior of 
dwarf galaxies and globular clusters. Again we find that Ret 
2 follows the galactic relation. Moreover, assuming dynamic 
equilibrium and negligible contamination from binary stars. 
Ret 2 has among the highest dynamical mass-to-lig ht ratios 
of any known object. The crude mass estimator off 


et al. (2009a i implies that the dynamical mass enclosed within 
Ret 2’s projected halfiight radius is M(Rh) ~ 5RhO’y,„/(2G) = 
2.4/Qg X 10^ Mq, and the associated mass-to-light ratio is 
« 2M(Rii)/Ly = 462/J 57 in solar units. 

Tablefflsummarizes the observed properties of Ret 2, based 
on the pnotometric results of K15 and the spectroscopic re¬ 
sults presented in this work. Where photometric results from 
DES15 differ from those of K15 (e.g., for absolute magni¬ 
tude), we list the DES15 results as well. 


5. SUMMARY AND DISCUSSION 

We have presented results from our initial stellar- 
spectroscopic observations of Ret 2. Integrating for two hours 
in below-average observing conditions, Magellan/M2FS has 


^ Many popular dynamical mass estimators have Mjyn oc /G, where 

the constant of proportionality is typically between ~ 2 —6 (e.g., |Richstone| 
|& Tremaine|1986[|Wolf et al.|2Q10[]Amorisco & Evans|2Qll| . 
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TABLE 4 

Summary of Reticulum 2’s observed photometric and spectroscopic properties 


Quantity 

value 

reference 

R.A. at center 

Q;j2ooo =03:35:42 

K 153 

Dec. at center 

<5j20oo = -54:02:57 

KI5 

Galactic longitude 

/ = 266.2958 deg 

KI5 

Galactic latitude 

^7 =-49.7357 deg 

KI5 

distance modulus 

m-M= 17.4 ±0.2 

KI5 

distance from Sun 

D ~ 30 kpc 

KI5 

absolute magnitude 

Mv =-2.7±0.1 (-3.6±0.1) 

K15 (DES15) 

exponential scale length 

Re = 3.37^'^3 arcmin 

KI5 

ellipticity 

e=l-(Ufl) = 0.59!jJ;Sj 

KI5 

position angle 

PA= 71 ± 1 deg 

KI5 

projected halflight radius 

Rh = 3.64 !^^2 arcmin 

KI5 

projected halflight radius 

«h = 32!ljpc 

K15 

systemic line-of-sight velocity 

Vios =64.8f] |!, km s“' 

this work 

systemic line-of-sight velocity 

v’los = -90.9 km s“^ 

this work 

internal velocity dispersion 

^^^ 0 ,= 3.61)!;’km s-' 

this work 

velocity gradient 

^^'los “ 2 ^ ^ * arcmin ^ 

this work 

PA of velocity gradient 

».l„s=-™1f deg 

this work 

mean metallicity 

([Fe/H]) =-2.67/!!;!) dex 

this work 

metallicity dispersion 

'^[Fe/Hl =0.50!2;l3 

this work 

metallicity gradient 

^[Fe/H] = dex arcmin”^ 

this work 

mass enclosed within 

M(«h)= 2.4/J;! X 10= M0 

this work 

mass-to-light raio 

T = 462/j“Mo/L0 

this work 


t Unless otherwise noted, K15 and DES15 report similar values. 


notes 


semi-major axis 


/?h 1.68/?^- 


solar rest frame 
Galactic rest frame, given solar motion measured by 


Schonrich et al. 


2010 


similar to median metallicity error 

M(/?h) ~ /(2G); assumes equilibrium, negligible binary stars 

T ~ 2M(/?h)/Lv; assumes equilibrium, negligible binary stars 


delivered spectra sufficient for estimating vios, Tgff, logg 
and [Fe/H] for 38 stars, reaching a limiting magnitude of 
r < 20.3. Based on a combination of photometrically- and 
spectroscopically-derived quantities, we have confirmed 18 
likely members of Ret 2. From the member sample we es¬ 
timate a velocity dispersion of 3.6'^'^ km s“* about a mean of 
64.8;!;[ g km s“^ in the solar rest frame, and a metallicity dis¬ 
persion of 0.50 ;!;q} 3 dex about a mean of-2.67;|;Q34 dex. While 
Ret 2’s size and luminosity are similar to those of globular 
clusters as well as ultrafaint dwarf galaxies, our estimates for 
mean metallicity and velocity dispersion place Ret 2 on scal¬ 
ing relations that are followed by dwarf galaxies and not by 
globular clusters. On this basis, we conclude that Ret 2 is a 
bona fide galaxy. 


Furthermore, under assumptions of dynamic equilibrium 
and negligible contamination from binary stars. Ret 2’s com¬ 
bination of size, velocity dispersion and luminosity implies 
a dynamical mass-to-light ratio of ^ A&I^^Mq/Lq, among 
the largest inferred for any known object. However, the va¬ 
lidity of these assumptions remains unclear. Perhaps under¬ 
mining the assumption of dynamic equilibrium is Ret 2’s sig¬ 
nificantly flattened morphology (K15, DBS 15), which—given 
Ret 2’s Galactocentric distance of ^ 30 kpc—may signal on¬ 
going tidal disruption. On the other hand, our velocity sample 
does not show a significant velocity gradient that mi ght reflect 
the ordered motions associated with tidal streaming (|Piatek &| 
Pryor|1995] Read et al.|2006| l. The mean velocity that we es¬ 
timate for Ret 2 becomes vios -90.9 km s“* in the Galactic 
rest frame, implying that Ret 2 is currently approaching peri- 
center. We expect that deeper photometric and deeper/wider 
spectroscopic followup will be required in order to determine 
the extent to which Ret 2 is undergoing a disruptive tidal en¬ 
counter with the Milky Way. 

Finally, even if Ret 2 is resilient to tides, its velocity dis¬ 
persion reflects some combination of its internal gravitational 
potential and the intrinsic vel ocity variability of binary stars. 
[McConnachie & C6te| (|2010|l have demonstrated that binary 


motions alone can contribute velocity dispersions of up to 
a few km s“'. Therefore we cannot exclude the possibility 
that Ret 2’s velocity dispersion receives a significant contri¬ 
bution from binary stars, in which case the dynamical mass 
and mass-to-light ratios that we have reported would be sys¬ 
tematically overestimated. More epochs of spectroscopic ob¬ 
servations of Ret 2—and other ultrafaint dwarfs with small 
velocity dispersions—are required in order to gauge the mag¬ 
nitude of this effect. 
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